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METHOD AND APP^ 



lGING THEOUGH 



SCATTERING OR OBSTRUCTING MEDIA 



Field of the Invention 

The invention relates to the field of ^..optical imaging. Specifically the 
invention relates to methods and apparatus for determining the optical 
temporal response of a medium to a short optical pulse excitation. 

Background of the Invention 

Pubhcations and other reference materials referred to herein, including 
references cited therein, are incorporated herein by reference in their 
entirety and are numerically referenced in the following text and 
respectively grouped in the appended Bibhography, which immediately 
precedes the claims. 

There are many appUcations that wotdd benefit from a noninvasive optical 
means of viewing objects embedded within turbid media or otherwise 
obstructed by other objects. In military and security appHcations the ability 
to see through obstructions or turbidity is of great importance for example 
for detecting the presence of concealed weapons under clothing or viewing 
enemy forces concealed by fog or cloud cover. 
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In medical diagnostics the most common methods for diagnosing 
abnormahties of internal organs are based on circumstantial evidence. For 
example, high cholesterol levels in the blood might be circumstantial 
evidence for vascular diseases. Obviously, the physician can use endoscopes 
to obtain a much better perspective of the extent of the disease. The ability 
to "see" the organ, i.e., to measure its optical properties, results in a wealth 
of clinical information. The reason for this is that most of the biological 
molecules absorb light in the optical regime of the spectrum. Therefore, it is 
relatively simple to distinguish between different biological tissues with 
optical devices. However, endoscopy and similar procedures are invasive 
methods and therefore cannot be used in many cases or. as frequently as the 
physician would like. 

For these reasons, noninvasive methods were developed to diagnose internal 
organs. However there exist few, if any rehable noninvasive techniques to 
measure the optical properties of these internal organs. As a consequence, 
most of the ubiquitous methods rely on circumstantial evidence for 
diagnosing the internal abnormalities. For example, in x-ray mammography 
most of the tumor findings are due to the presence of micro-caldfications, 
which gives only circumstantial evidence of the tumor's presence. 

The inability to see through biological tissue is only a specific example of a 
broad family of related problems pointing to the apparent futihty of using 
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optical radiation for direct imaging thxougli a turbid or obstructing medium. 
This problem is ubiquitous in everyday life. Fog, clouds, t\irbid water, and 
milk are just a few examples from everyday encounters with such optically 
tuxbid media. In particular, biological media are hi^y diffusive for visible 
and infrared light. As a result, an optical image is severely deteriorated, 
even after passing through a few millimeters of tissue. 

Despite this severe shortcoming, it is clear that a reliable non-invasive 
optical technique would be a very useful clinical tool in diagnosing in-vivo 
internal organs. Moreover, the use of light as a diagnostic tool has many 
advantages over other non-invasive diagnostic techniques. Unlike x-rays, 
light is not oncogenic, it has the potential of giving a much higher resolution 
than ultrasoimd imaging, and the required equipment is potentially much 
cheaper than MRI. These benefits and technological progress in optical 
equipment have encouraged researchers to develop optical diagnostic -tools 
for turbid biological systems. 

Although many methods based on coherence, polarization, diffusion 
properties etc. have been developed [1], two techniques appear to be the 
most promising: the time-domain and frequency-domain techniques. [Other 
techniques, such as optical coherence tomography, which utilizes short- 
coherence-length Hght, can be regarded as equivalent to the short p\iLse 
illtunination method] 
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The most direct and one of the most popular approaches to seeing through a 
turbid medium is the first-light or time-domain technique. In this method 
the first arriving light, which contains the undistorted image, is separated 
from the scattered light by a very fast time-gating camera [2-6], In ref. 3, 
the authors suggest that a resolution of a few millimeters is achievable in a 
system with a temporal resolution of about lOps. This method demands 
complex and very expensive equipment. 

While the time-domain technique gives superfluous information, and the 
requirement of the measurement equipment is to distinguish between the 
informative (i»e., ballistic) portion of the signal and the noise (i.e. diffusive 
portion), the frequency-domain techniques usually resvilts in less 
information. In the frequency-domain approach [1,7-9] the laser light is RF 
modulated (MPIz-GHz), and the measin-ement of the amplitude and phase of 
the modulated diffusive Hght take place at the boundaries of the medium. 
Then, a numerical reconstruction takes place in which the hght distribution 
in the medium is calculated from the collected data. 

The frequency-domain technique is considered to be simpler, more reliable 
and more economic than tiie time-domain method, since it needs much 
simpler equipment. Moreover, unlike the time-domain method, the detected 
signal in this technique can be considerably stronger since it includes the 
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entire signal energy, not merely. the ballistic or quasi-ballistic portion, which 
can be quite small in a Mghly turbid medium, 

In principle, the two methods are eqtdvalent since the latter one can be 
regarded as the Foiirier coiinterpart of the former. However, except for 
small technical details, the main difference between the two is the spectral 
range. In the time-domain techniques (for medical purposes) pulses of about 
a few picoseconds are used, which is equivalent to a spectral range of more 
than ipOGHz. However it is difficult to modulate light at frequencies 
greater than lOGHz, Therefore the concept is different: instead of collecting 
the data in time, in the frequency- domain method the data is collected in 
space (on the medium's boundaries). Such a method requires complicated 
numerical algorithms, is sensitive to boundary conditions, and its spatial 
resolution is hmited by the medium's diffusion length. 

It is therefore a pvirpose of the present invention to provide a method for 
determining the optical temporal response of a medium to a short optical 
pulse excitation, at least a portion of which is indicative of the position and 
shape of an object or objects embedded in an optically turbid meditun and/or 
obstructed by other objects that are at least partially transparent. 



It is another purpose of the present invention to provide a method for 
determining the optical temporal response of a medium to a short optical 
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pulse excitation, at least a portion of which is indicative of the position and 
shape of an object or objects embedded in an optically txxrbid medium and/or 
obstructed by other objects, that can be carried out without using short 
pulses and fast detectors. 

It is a further purpose of the present invention to provide a method for 
determining the optical temporal response of a medium to a short optical 
pulse excitation, at least a portion of which is indicative of the position and 
shape of an object or objects embedded in an optically turbid medium and/or 
obstructed by other objects, that can be carried out with relatively 
inexpensive equipment. 

Further purposes and advantages of this invention will appear as the 
description proceeds. 

Sutmnarv of the Invention 

The present invention relates to providing a method, which wiU be called 
herein Spectral Ballistic Imaging (SPEBI), of computationally deriving the 
ballistic photon (or "first ]igh^') image after characterizing the optical 
spectral transfer function of the medium in a suitably wide spectral band 
and with suitable spectral resolution. Let /^(/) describes a short pulse 

which, in the prior art method, would be sent through the medium in a 
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^ballistic' imaging system, and f^^it) is the output which, would be obtained 

and would contain the information on the ballistic and 'first light' response, 
as well as the light which would arrive later due to scattering. In the 
present invention, the Fovirier transform jE^(a?)of the short input pulse is 

sent instead. The output spectrum (o) is measured, giving the medium's 
transfer function By determining this spectral transmission response 
of the medium with a suitable resolution, an inverse Fourier transform 
operation will allow the determination of the short p\ilse response /o„,(0 In 

this way the *fate' of the ballistic photons can be determined and a high- 
resolution image can be achieved. This invention allows for trans- 
illumination imaging through a scattering volume, or reflection-illumination 
imaging through discrete scattering layers. 

In a first aspect, the present invention is directed towards a method for 
determining the optical temporal response of a medixim to a short optical 
pvilse excitation. The method comprises sending hght comprising spectral 
frequencies, which make up the Fourier transform of the short pulse to 
be emulated through the medium. The relative amplitude and phase 
change of each of the spectral components hght exiting the medium is 
determined with respect to that of the illuminating light source and the 
spectral response of the medium is obtained from the relative amphtude 
and phase change. An inverse Fourier transform is then 
computationally performed on the spectral response to obtain the 



wo 2004/008116 



-8- 



PCT/IL2003/000582 



temporal response of the medium to the emulated short pulse. The hght 
can be either CW (Continuous Wave) or modulated. The method can be 
carried out using light that comprises only part of or substantially all of 
the spectral frequencies which make up the Fourier transform of the 
short pulse to be emulated. The method can be carried out by using light 
which has a carrier frequency that is scanned over time; and the relative 
amphtude and phase change are determined for each carrier frequency. 
The method can also be carried out using hght consisting of a broad 
spectral bandwidth and each of the spectral components of the output of 
the hght exiting the medium is detected. 

The hght can be detected either at a point, along a line, or over a two- 
dimensional area. The medium can comprise discrete layers. The detected 
hght can be either transmitted through the medium, reflected back from the 
medium, or exit the medium at any angle with respect to the illumination 
path. There can be at least one object embedded within the medium which 
can be detected and imaged. The at least one object can comprise a 
cancerous growth or other undesirable physiological formation, and the 
medium comprises a segment of the human body..' 

The *first light' response, or other specific segment of the temporal 
response, which is indicative of the position and shape of the object or 
objects embedded in an optically turbid medium and/or obstructed by 
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other objects, which are at least partially transparent, cab be determined 
from the optical temporal response. 

In another aspect, the present invention is directed towards a system for 
determining the optical temporal response of a medium to a short optical 
pulse excitation. The system is comprised of a timable laser, an RF 
oscillator, modulatiag means, detecting means, an optically scattering 
medium, electronic processing means, and optional optical elements 
means. 

All the above and other characteristics and advantages of the invention will 
be further imderstood through the following iQustrative and non-hmitative 
description of preferred embodiments thereof with reference to the 
appended drawings. 

Brief Description of the Drawings 

- Fig. 1 schematically shows the SPEBI system in a transmission 
configuration; 

- Fig. 2 schematically shows the SPEBI system iu a reflection 
configuration; 

- Fig. 3 schematically shows the modulation and frequency scanning 
scheme for a preferred embodiment of the SPEBI system; 
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- Fig. 4A schematically shows the experimental setup for a SPEBI 
experiment to detect the balhstic hght response through diffuse glass 
and a hidden object; and 

- Pigs. 4B and 4C show the results of the SPEBI experiment of Fig. 4A. 

Detailed Description of Preferred Embodiments 

In this specification the terms ''turbid medium'', "scattering medium", 
"optically scattering medium", "optically tin-bid medixmi", and other terms 
having a similar meaning are used interchangeably. 

Fig. 1 describes a transmission system used for seeing through a turbid 
optical medium to acquire an image of an object obscured by the medium. 
In the preferred embodiment shown in Fig. 1, a CW frequency-tunable laser 
10 is modulated at an RF frequency Q through the use of a RF modulation 
source 17, which modulates the laser with the electronic sinusoidal signal 
18. In the weak modulation regime, the modulated hght field 11 is 



(/, <0j ) = e""^' {l + a cos Qt) = 



2 ^ 



[lA] 



where coj is the optical carrier frequency and a«l (for simphcity, unit 

amphtude at the carrier frequency is assumed). It should be noted that the 
weak modulation regime is neither an essential nor a simplifying 
assumption; similar conclusions are obtained for deep modulation, where 
expression (1) takes a simpler form 
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E^(t,<Oj)^e'^^'(e'^+e-"'). [IB] 

Going back to the weak modulation case, the intensity is 

//^(/) = l + 2acosQ/. [2] 

where higher orders of the small term a were neglected. 

A portion of this beam is directed through the turbid medium 12, which 
includes an embedded object 13. The output field 14 is: 



2 



[3] 



The added phases cpj-, ^/*^and ^j^"^ for the three spectral components at cOy , 

fljy + Q and a)j-Q respectively, are due to the effective optical path length 

difference at these three frequencies. These path length differences are 
mainly due to material dispersion and frequency-dependent scattering 
paths. The factor A((d.j) represents frequency-dependent attenuation and is 

real-valued. For simphcity, it is assumed for now that this factor is constant 
over the frequency range a?y ±Q, The phases as well as the amplitude can 

also change due to the electronic system, cables, detector, etc., however 
these factors are independent of the scattering medium and embedded 
objects, so that they are determined and accounted for. Therefore, the 
output intensity at 14 which is detected by detection system 15 is 
s ^2(a)y').[l + 2flccos(Q^ + A9((i)J)], [4] 
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where 



(+) 



2 



and c s cos (pj 



2 



[5] 



J 



Electronic processor 19 compares EF signals I^^ (Eq.4, 16 in the figure) to 

(Eq.2, 18 in the figure), so that the ampHtude and phase values A(oi>j) 

and ^pia>j) (20 in the figure) can be determined. 19 can be a phase detector, 

network analyzer, lock-in amplifier, or other electronic processor known in 
the art. 

The laser fi-equency is tuned in steps of Sco over a frequency range A(d and 
the values of A{a)j) and Ag){cDj) are acquired for every frequency coj . 

Fig. 3 schematically shows the modulation and frequency scaiming scheme 
for a preferred embodiment of the SPEBI system. In the figure the top line 
30 shows the phases of the jth step and the bottom line 31 the phases of the 
0' +l)th step. The arrows 32 show that cp^^^ = q)^~^j for every step y. 

If the frequency steps Sa> fulfill the relation shown in Fig. 3, i.e. 



then, the spectral response of the medium, within the range Ao) , can be 
determined unambiguously: 



Sco = 2Q 



[6] 
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H(ci)j)-A((Oj)e' 




[8] 



m=0 

is the accumulated phase, and (D^ is the lowest laser frequency to be used. 
If the frequency steps are larger, i.e. 5© > 2fl , Eq.(8) should be interpolated; 



In the above derivation, it is assumed that the ampUtudes A{cOj) are 

constant over the frequency range (Dj±Q. This is not a necessary 

requirement. If they differ, then the output intensity (Bq. 4) will have a dc 
value and a peak-to-peak ac value, which is dependent on the ampUtudes, so 
that the ampUtude values can be determined from a measurement of the dc 
and peak-to peak values. 

Referring to Eq. 4, it is to be noted that the phase variable c , which affects 
the modulation depth of the modulated light at the output, is also dependent 
upon the added phases of the spectral components. Therefore, the phases 
(pifDj j can also be determined by measuring the modulation depth of the 

output signal 
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A computational inverse Fourier transform of H{a)) will give the temporal 
response of the medium to an input pulse whose spectral bandwidth is Afi> . 
The 'first-light' signal can be extracted from this information by 
computationally filtering out the bxolk of the temporal response and 
retaining the 'first-hght' response (in analogy to the task of the fast shutter 
in a real short-pulse measurement). Note that the computational work, 
which may form the hulk of the acquisition time, can be done 'off-line*, either 
, simultaneously, or after the spectral data is collected. 

Other holographic methods [10] are in some sense similar to aspects of the 
current invention, since in both methods the measurements are done in the 
spectral domain, and both amplitude and phase are measured. However, 
there are two crucial differences between the two techniques: 1) The 
holography methods are based on interference of two beams that pass 
through two totally different trajectories. This type ' of interferometric 
measurement is highly sensitive to system noise, and difficult to implement 
in non-laboratory systems. In the present invention, on the other hand, the 
two beams copropagate along common paths, and therefore the 
measurement is extremely robust; 2) Moreover, since in both methods the 
phase and ampHtude must be measured, the methods of Arons et al [10] 
consists of two scans: wavelength and time delay, while the present 
invention requires only a wavelength scan, since the measinrements of 
ampUtude and phase are done simultaneously. 
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Fig. 2 shows schematically a preferred embodiment for acquiring the image 
of an object obsctired by an absorbing layer in a reflection-illumination 
configuration. Elements 10, 11, 15, 16, and 17 to 20 are identical to those of 
Fig. 1. In Fig. 2, the laser output 11 is directed towards a target, which 
consists of multiple scattering layers, two are depicted in the figure as 22 
and 23. Layer 23 is also the object, hidden behind the front scattering layer 
22. A portion of the Ught 21 transverses the first layer and is reflected back 
by the object at small enough angles with respect to the optical axis that it 
is collected by the detection system. Much of the light is scattered in a wide 
angle by all of the layers, shown as 26. A beam sphtter 24 deflects the. 
returning light 25 into detector 15. 

The laser modulation and signal processing is identical to that of Fig. 1. As 
in the transmission case, the temporal response is computed, and the 'first 
Hght' response is determined, so that the image information of the hidden 
object can be extracted. 

As in time-gated imaging, as well as other imaging techniques, it is difficult 
to extract information on objects deeply embedded within a scattering 
voliime when using the reflection configuration. This is due to the 
accumulated backscattering of hght from the scattering volume in front of 
the object, which is usually much stronger than the ballistic backreflection 
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from the desired object. Therefore, SPEBI, hke other first-hght techniques, 
works best with discrete scattering layers, as shown in Pig. 2, or in 
situations where the objects are not deeply embedded within the scattering 
volume. The balhstic hght, which is reflected from each layer would appear 
temporally separated, so that the information on the object layer can be 
extracted. 

Acquisition of a complete two-dimensional (2D) image of the desired object 
requires the computation of the ballistic hght components within the entire 
3D volume. This can be done by scanning the input light and detector 
throughout the 2D input and output planes, and performing the SPEBI 
measurements for each coordinate. Alternatively, the .entire input plane of 
the scattering medium can be illuminated simultaneously, while performing 
parallel measurements using a detector array at the output plane. The hght 
sources and/or detectors may include lenses, optical fiber, mirrors or other 
optical components for hght expansion, collection, focusing, steering, etc. as 
is weU-known in the art. 

The above discussion leads to the conclusion that the main parameters to be 
determined to enable the SPEBI measurements to be carried out are the 
desired spectral resolution 5(o, laser modulation frequency Q and the 
required tuning range A^. These parameters and the method used to 
determine them will now be discussed. 
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The tuning range is optimally the spectral bandwidth of .the input pulse 
width , which is to be emulated with the txmable CW laser. A good rule of 

thumb is that for a dense optical medium of length Z, the pidse width should 
be significantly smaller than the 'flight time' of a ballistic photon - nil c 
through the medium, so that r^^Q.Olnljc , where n is the effective index of 
refraction. Assuming a Gaussian pulse, so that Afi? = 4/r^, gives 
Lay = 400c/ nl . The required spectral resolution ScD is dependent upon the 
characteristic diffusion time of the medium. In order to obtain full 
reconstruction the spectral resolution of the measurement should be better 
than dco ^Itt I Tj^ where is the diffusion time, which can be 

approximated by = / IDB'ifi^'nl /c.Dis the diffusion constant and 
III/ is the effective scattering coefficient (the reciprocal of the random walk 
length). As an example, for ordinary human tissues 3^ = 10ns .' 

Assuming 100 resolution points, or Aa)/So)>10^j gives &y = 4c/«/ and 
Q = 2c/ nl . Assuming, for example, a length of 10 cm and w = 1 .5 , leads to the 
following required operating parameters for a tunable laser with a center 
wavelength of 1530nm: a wavelength tuning range of l.Onm, timing 
resolution of 10pm, and modulation frequency of 640 Mhz. This will 
emulate the performance of a 5ps laser pulse. 
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Tunable lasers having these parameters are commercially available. For 
example, the TUNICS-PRI-1530 supplied by NetTest Photonics fulfills the 
above requirements. This laser has the following specifications: wavelength 
txinability range 1480-1580 nm, wavelength resolution =lpm, and power 
modxilation firequency Q/2;r up to 1 GHz. The tuning range of lOOnm 
means that a minimum pulse length of approx. 50 fs can be emulated. The 
1pm resolution capabihty means that the lowest 5o/2it = 128Mfe and that 
up to 10^ resolution points can be measvired. Skilled persons wiU realize 
that the number of resolution points can be lowered to a Hmit where 

An important advantage of the timable laser is the ease with which its 
power can be amplified. Since it operates at the optical communication 
wavelength band and is fiber-coupled, its power can be easily (and relatively 
cheaply) amplified from 3 mW to over 500 mW using an Erbium-doped fiber 
amplifier (EDFA). For medical and other applications, it may be 
advantageous to employ wavelengths in the 800nm-900nm band. Other 
types of timable laser sources, e.g. Ti:Sapp,, diode lasers, that are known in 
the art can be used in SPEBI. 

Alternatively, a broadband hght source of bandwidth b.co can be used to 
illuminate the sample; while, at the output, the detector system includes a 
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means of spectrally separating the light components and meas\rring the 
amplitude and phase of each spectral component. 

The actual number of resolution points that will be required in a given type 
of mediiim may be lower, and in certain circumstances significantly lower 
than the result of the above model. For example, if a priori information 
regarding the medium^s optical characteristics is known, then the number of 
resolution points can be decreased. In many circumstances, interpolation of 
the missing data woxald be sufficient. In certain applications, a data bank of 
spectral data, e.g. data on breast tissue with and without cancerous 
growths, could be used in order to efficiently categorize the data of a 
particular patient, so that a hig^-resolution measurement would not be 
reqviired. 

In Fig. 4A is schematically shown the SPEBI set-up 40 for an experiment 
designed to measure the time response . of opal diffuse glass 12, behind 
which was place a small absorbing object 13. The glass was illuminated 
with light from a timable laser 10 of the type described above. The effective 
input pulse length was 1.5ps, the modxilation frequency was IGhz, and the 
thickness of the glass was 3mm. Also shown in Fig. 4A are RF oscillator 41, 
modulator 42, optional optical elements 48, detector 15, and electronic 
processor 19 which provides the amphtude and phase data. 
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Fig, 4B shows the reconstructed time response of the medium, with and 
without the absorber, to an input pulse of 1.5ps shown at t = 0. According to 
the method of the invention, it is not this pulse, but its Fourier transform 
that is transmitted through the diffusing medium and the object. The 
inverse Fourier transform of the recorded data is then calculated and it is 
this that is known as the reconstructed optical response of the diffusing 
medium and the absorbing object which appears between 0.5 and 1.1 x lO-i^ 
seconds after the initial pulse. 

Fig. 4C is an enlargement of the features shown between 0,5 and 1,1 x lO-^^ 
seconds in Fig, 4B. In Fig. 4C, curve a is the reconstructed optical response 
of the medium alone and curve b that of the medium with the absorber. 
Curve b has been magnified by a factor of 7. In the reconstructed optical 
responses, pulse expansion and deformation due to the scattering medium is 
clearly seen. In addition, it can be seen how the absorber has blocked the 
first-arriving hght. This experiment demonstrates that a temporal 
resolution of less than -^Ips is easily achievable; that it is possible to 
separate the baUistic and quasi ballistic ("snake photons'') portion of the 
signal from the diffusive component; that an obstructed object can be 
detected; and, perhaps most importantly, that a relatively noiseless signal 
can be obtained with this technique. 
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Although embodiments of the invention have been described by way of 
illustration, it will be understood that the invention may be carried out with 
many variations, modifications, and adaptations, without departing from its 
spirit or exceeding the scope of the claims. 
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